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Abstract: We show how to correctly account for scalar accretion onto black holes 
in scalar field models of dark energy by a consistent expansion in terms of a slow roll 
parameter. We find an exact solution for the scalar field within our Hubble volume 
which is regular on both black hole and cosmological event horizons, and compute 
the back reaction of the scalar on the black hole, calculating the resulting expansion 
of the black hole. Our results are independent of the relative size of black hole and 
cosmological event horizons. We comment on the implications for more general black 
hole accretion, and the no hair theorems. 
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Black holes are one of the most fascinating objects in classical general relativity. 
They represent the endpoint of gravitational collapse of a large mass object, indepen- 
dent of initial conditions. The prototypical black hole, the Schwarzschild solution, 
was first presented barely a year after Einstein put forward his new theory of gravity, 
yet it took half a century before relativists were confident of the interpretation of this 
solution, and began to construct a rigorous set of theorems describing the properties 
of black holes. One of the more poetically labelled set - the "no hair" theorems []J 
has been perhaps the most contentious. Originally demonstrated for static vacuum 
black holes, p|, [|, the no hair theorems have been refined and extended to a wide 
range of interacting particle and gravitational theories, @, |[ |J . Nonetheless, there 
are many known "violations" of the no hair theorems, from the dressing of black 
holes by scalar or gauge condensates |7|], to literal hair, in the form of topological 
defects which extend out from the black hole to infinity, ||. 

The issue of the most basic hair however, scalar hair, is thought to be understood, 
yet is perhaps the most perplexing. For a wide range of potentials, it is easy to show 
that a static black hole cannot have a nontrivial scalar field on its horizon (unless 
the spacetime is asymptotically anti-de Sitter, |§, or the scalar potential has been 



specifically tailored, [10]). However, within cosmology, scalar fields are widely used, 



not only for inflation, but as an expedient model of quintessence, [11]], or late time 



acceleration (see JT^j for a review of dark energy models). Indeed, the standard 



exponential scalar potential - widely used to give power law acceleration, was shown 
to be in conflict with spherically symmetric asymptotically flat or de Sitter static 
black hole spacetimes, ||13|| . How then can these cosmological scalars interact with a 
black hole horizon? Generally, the scalar field is varying on a very slow time scale 
relative to the light crossing time of the black hole, thus we might expect that it has 
little effect on the horizon itself, yet, like the inexorable flow of time, the scalar must 
evolve cosmologically, and it seems contradictory that it is pinned to a single value 
for all time at the black hole. 

In fact, the resolution of the conflict is quite straightforward: The no hair the- 
orems were explored in the context of static black holes, and quite clearly a rolling 
cosmological scalar can never be static. Correspondingly, the presence of a black 
hole in a cosmology implies that the cosmology cannot be spatially homogeneous, 
and therefore does not have the canonical FRW form. While we can use cosmolog- 
ical perturbation theory to describe the far field effect of the black hole, we cannot 
describe the full spatial geometry including the event horizon within the context of 
perturbation theory - and it is the effect of the event horizon which is critical in 
determining the behaviour of the scalar. We therefore need a time- dependent black 
hole solution. 

Exact black hole solutions are surprisingly thin on the ground, given how ubiq- 
uitous black hole are in nature, and time dependent ones even more so. The first 
dynamical black hole 'solution' was postulated by McVittie, |T5| , who modified the 
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Schwarzschild solution to include an FRW-like cosmology, with the black hole horizon 
at fixed comoving radius and a fluid energy density dependent only on time. However, 



as pointed out in fL6|| , one expects the black hole to break the spatial homogeneity 
of the cosmology, and thus (except in the special case of a de Sitter universe) this 
unrealistic set-up leads to a singularity on the would be horizon. The Vaidya solu- 



tion, [n], also time dependent, is somewhat special, as the mass of the black hole 
now becomes dependent on a null coordinate, requiring a null matter source. In 
cosmology however, we expect a black hole to accrete general stress-energy, and as 
already noted, solutions in the literature with nontrivial scalar fields have very spe- 
cialized potentials, with fixed mass black holes. Thus a direct approach to finding 
an exact solution seems to have been unsuccessful in that solutions do not account 
for an evolution of the black hole horizon area due to scalar accretion. 

An alternate approach in the literature is to find a probe solution for the scalar 



field, i.e. one in which the scalar field evolves on a Schwarzschild background, [[17, 



18fl . Here, near the horizon, the scalar field must be a function of the advanced 
null coordinate, v ~ e t+r , which parametrises the black hole future event horizon. 
One can then estimate the back-reaction on the event horizon using the energy 
momentum of this approximate solution. A problem with this approach however, is 
that the coordinates being used are local static Schwarzschild coordinates, which do 
not correspond to the cosmological flow of time at large distances and therefore give 
no guarantee that any 'solution' will be well behaved at cosmological horizons. One 
can also worry that estimates using local intuitive notions of energy, rather than an 
analysis of the Einstein equations with a scalar source, might be misleading. 

In order to be confident that probe calculations give a good estimate of scalar 
accretion onto black holes, we require a time- dependent scalar field with a time 
dependent black hole. Here, we present a resolution of this problem, by expanding 
the equations of motion for the geometry and the scalar field order by order in a "slow 
roll" parameter. We show how to put a time dependent scalar on the black hole, 
and how this extends to the cosmological solution at large distances. Our solution is 
fully nonperturbative in the radial direction, thus is valid independent of the relative 
sizes of the black hole and cosmological horizons. We also present the back reaction 
of the scalar field on the black hole, and calculate the expansion rate of the black 
hole due to scalar accretion. An advantage of our method is that we can identify the 
event horizons accurately as null surfaces, therefore the usual ambiguity of apparent 
vs. event horizons does not occur. 

To motivate our approach, consider a universe in which acceleration is driven by 
a scalar field with an exponential potential W((f>) = M 4 e _/3</ \ Writing the solution 
in conformal time gives 



ds u 
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\[drj — dp)(dr/ + dp) — p 2 dVt\ ] \ 
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= 0o-^^yln^ (2) 

where e = (3 2 /2k = 2W(s) 2 * s none other than the slow roll parameter for the 
potential W, and represents the power of the acceleration. Clearly, if we assume that 
£< 1, [Vo/v] 1 ^ ~ 1 — 2e In 1 77/770 1, and hence this metric is perturbatively close to 
the de Sitter (dS) metric, over a Hubble time interval. From (|j) , <f> = </> + O (e/k) ' , 
and thus our cosmology is roughly dS, with a small correction of order 0(e 1 ^ 2 ) for 
the scalar, and 0(e) for the geometry. 

Now consider adding a black hole to this cosmological rolling scalar solution. 
The geometry will not now have constant curvature spatial slices, but we do expect 
an SO (3) symmetry, corresponding to the spherical symmetry of the black hole. We 
also expect both time and radial dependence in the metric. The general metric may 



be written in the form |L9], ^0 



ds 2 = Ae 2u B- 1/2 dUdV - Bdn 2 u . (3) 

This form of the metric elucidates both the symmetry of the spacetime (SO (3)), as 
well as the remaining gauge freedom (the conformal group in the U, V directions). 
By rewriting the coordinates in lightcone form, it will be clear how to deal with the 
horizons present in the anticipated solution (the cosmological and black hole event 
horizons), as well as how to change gauge to analytically extend across these horizons. 
In particular, the null coordinates allow us to identify the actual event horizons of 
the solution, as opposed to apparent horizons, as an horizon is of course always a 
null surface, defined by U or V = constant. 

The coupled black hole scalar equations we have to solve are 

B tUV = 2 {kW{<P)B 1 / 2 - B- 1 ' 2 ) e 2v (4) 

ujuv = \ {kWWB-W + B-*' 2 ) e 2 » - ^My (5) 

4>,uv = -W^B-^e 2 " - ± {B^ y + B iV (f> }U ) (6) 

By V = IVyBy - K B<f?y (7) 

B.jju = 2u tU B tV - kB4? v (8) 

For a constant scalar field (as one might naively suppose from the no hair theorems) 
the Einstein equations have Schwarzschild de Sitter (SdS) as a general solution, [|Kf|. 
To briefly review, if is constant, (|7|) and (|^) can be integrated directly to show that 
B must be a function of [-F(V) + G(U)], where F and G are arbitrary integration 
functions, and e 2u = F'(V)G'(U)B'(F + G). Equation (§ then gives 

B' = -kW B 3/2 - AB l/2 + 8GM = -4v / SA^(v A B) (9) 

o 
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where SGM is an integration constant (suggestively labeled!), and we see that the 
function N defined above is none other than the SdS potential, with kWq/3 = H 2 , 
the Hubble parameter. One can then change coordinates to r = I? 1 / 2 , t = 2{G{U) — 
F(V)), to obtain 

F"G' /?' dr 2 

ds 2 = „° AdUdV - BdQl = N(r)dt 2 - -^r ~ r 2 dSll (10) 

B 1 / 2 u v ' N(r) u v ; 

i.e. the Schwarzschild de Sitter metric in static coordinates. 

Although the static coordinates are useful to identify horizons in the SdS metric, 
they are distinct from the "cosmological" coordinates, in which the metric would 
at large r asymptote a standard cosmological de Sitter spacetime; they also have 
the disadvantage that they are singular at both the black hole, and cosmological, 
event horizons. Writing r^ (r^ < < th < r c ) as the roots of N(r), and Ri = 
r.j/(l — 3H 2 r 2 ) = l/N'fri), we have a locally regular null coordinate system at each 
horizon given by 

V = R c exp [-2F/R C ] , U = R c exp [2G/R C ] CEH 

v = R h exp [-2F/R h ] , u = -R h exp [-2G/R h ] BH 

the latter being the standard Kruskal coordinates for the SdS black hole, and the 
former a set of Kruskals for the cosmological event horizon, provided we identify 
F = —(t + r*)/4, and G = (t — r*)/4, with r* being the usual tortoise coordinate. 

Now suppose that <p is not constant, but slowly rolling, then we can solve our 
spacetime and scalar field order by order in the slow roll parameter. The solution will 
be perturbative in e, but, crucially, nonperturbative in the radial coordinate across 
our Hubble volume. Writing (ft = (f> + 2e4>i/ (3 in (|) gives 

<f> 1>uv = m 2 B- l / 2 B'F'G' - %- {G'<j> iy + F'4> ltU ) (12) 

It is reasonably clear that the form 

(f) 1 = 2\(G(U)-F(V)) + <p(F + G) (13) 

will satisfy (J2J) provided [B<p']' = ?>H 2 B 1 I 2 B' ) hence 
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C 
R i \n\r-r i \ + ——\nr (14) 



with C an integration constant. 

For a nonsingular solution, the field must be regular in a locally regular coordi- 
nate system at both the black hole (r^), and cosmological (r c ), future event horizons, 
thus we get two constraints on A and C, solved by 

C = -H 2r jdpl±lA , A = (r r r ; } (15) 

r h + n r£ + r 2 
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Pulling this together we can write the <fi field in the static SdS coordinates: 

2e 



4 » + 7T A 



t — R c In \r — r c \ 



R h ln(r - r h ) 

Rn ( TT o C \ . . . r h r c In r 

--\ HrN + ^) Hr + n + rh) -V^5 

which is manifestly nonsingular at the horizons, and plotted in figure Q. 

SINGULARITY FUTURE INFINITY 




U+V=const 



t=const 



(16) 



Figure 1: Contours of constant <j> (thick blue) within our Hubble volume. For comparison, 
contours of constant t and r are shown in grey, and a constant cosmological time contour 
(U + V) shown in red. The horizons are labelled, as are the singularity and future infinity. 



It is useful to note the asymptotic forms of the scalar field. Near the cosmological 

event horizon 

■>~ 

(17) 



2e, 



^ + —2XR c \n(U + V) 



which gives the correct cosmological rolling scalar solution to order 0(e) as r% — > 0, 
(cf (H)). On the other hand, near the event horizon, 

0~0o + ^2Ai4hi?; (18) 

in agreement with the linearized solutions of [[17], [18]. Clearly, <p is regular from 
the black hole event horizon out to the cosmological event horizon, with regular 
derivatives, and thus its energy momentum is finite in this region. We can therefore 
compute the backreaction on the geometry to get a consistent solution to order 0(e). 
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Writing B = r 2 {l+e5i) and v = vsds+£$2i the procedure follows the integrability 
argument for the SdS solution: First integrate (0) and @ to obtain two expressions 
for 5 2 

2 ^2 = *i + r -^r + 9(U) +I v = Si- ^# + f(V) + Iu (19) 



where 



'v 



2R r N * y ' 2R C N 

W !*(*>, + £ -A) 



B v J N 2X r 2 

Next, consistency of ( |T9| ) determines the general form of 8\\ 



(20) 



HUV)(I V -Iu)- J ydV + J ^jrdU + h(r) (21) 



R C N 

with (|j) giving a second order ODE for h(r), solved by 

dr 



h(r) 



rN 2 R r 



hi 



[[l2H 2 r 2 <p + (rNy(Iu + I v )} 



(22) 



Regularity at both horizons then uniquely determines the various integration con- 
stants (or functions) and hence SB, 8v. 

The general expressions for SB and Su are rather cumbersome, but by focussing 
on the event horizons, it is straightforward to obtain results of most physical interest. 
For the cosmological event horizon, the {U, V} coordinate system is appropriate, with 
the CEH being at V = 0, and parametrized by U . Since r = const, along the horizons, 
(§) gives 

B = r 2 (1 - 8e\ 2 R 2 c In \U\) ~ r 2 \U\- 8£x2R " (23) 

Setting rh — > gives 4\ 2 Rl = 1, and hence B oc U~ 2e , in complete agreement with 
the cosmological horizon area of the pure rolling scalar solution, (|ip. 

Of more interest however is the accretion of scalar field onto the black hole. Here, 
using the Kruskal {u,v} system and (0), we get 

B = r 2 h (1 + 8e\ 2 R 2 h In v) ~ r 2 h v^ R l (24) 

We can compare this with an order of magnitude estimate based on naive physical 
notions of mass and energy flow, [17, 18|. The flow of energy into the black hole 



should be governed by the difference of the energy momentum tensor from being 
null, Tq — T r r , which is of order 2 ~ W' 2 /H 2 ~ eH 2 /k. Integrating this over 
the horizon gives k5M ~ r 2 eH 2 , or using the relation between horizon radius and 
mass: 5A ~ r\eH 2 5t. Of course, on the horizon, t is not a good coordinate, and we 
should instead use advanced time: t — > 2r h \iav, which gives 5 A ~ Ae(5\n.v)rf l /r 2 , in 
agreement with 
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For an astrophysical black hole, this accretion rate is glacially slow, and far 
outweighed by the local environment, in which the accretion disc far outweighs local 
interstellar matter, let alone this cosmologically coasting scalar. However, the fact 
that naive local estimates of the backreaction of accretion of scalar matter in this set- 
up are fully backed up by this radially nonperturbative analytic calculation means 
we should have confidence in these physically motivated techniques. 

The other feature we should emphasise about our results is precisely their non- 
perturbative nature in the radial coordinate. While the expansion is perturbative in 
time, our solution for the scalar, and hence the geometry, is exact across the space- 
time between the black hole and cosmological event horizons. Thus, we can apply 
these results to black holes whose event horizon constitutes a significant fraction of 
the Hubble volume. These black holes accrete at a similar rate to the increase in 
area of the cosmological horizon. 

We should point out that while we have a time dependent scalar field in a time 
dependent cosmological black hole background, this should not be viewed as a vio- 
lation of the no hair theorems. The solution corresponds to a rolling cosmological 
scalar field in which there is a black hole, and the rolling of the scalar adjusts to 
its presence. The scalar does roll in the vicinity of the black hole, although it lags 
behind the cosmological evolution, in the sense that the constant <p contours lie in 
front of the constant 77 contours in figure p]. However, the black hole does not have 
any scalar charge - there is no 1-parameter solution for in the black hole back- 
ground. The solution ( p!3|) is not the most general solution, there will be wave like 
fluctuations around this background, but we have not been able to find a family of 
solutions with spacelike dependence, which would be a signature of a scalar charge 
on the black hole. Thus, our results should be viewed as a way of reconciling the "no 
hair" intuition with more general time dependent situations. 

Finally, to our knowledge, this is the first analytic spacetime presented with a 
time dependent accreting black hole that results from consideration of a physically 
realistic matter system (rather than reverse engineering from a given metric). The 
results apply to a general potential, and only require that the scalar field is slowly 
rolling. As such, they represent a testing ground for investigation of black hole 
phenomena in the time dependent regime. 
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